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We investigate the enhancement in the D^D^tt^ final state with the mass M 

-0.3 



Dh; 3875.2 ± 0.7j:";g ± 0.8 MeV found recently by the Belle Collaboration in the B 

■ KD^D^TT^ decay and test the possibility that this is yet another manifestation of the 

^SJ . well-established resonance X(3872). We perform a combined Flatte analysis of the 
> 

^ ; data for the L>°D°7r° mode, and for the vr+vr- J/^J mode of the X(3872). Only if the 

^ ; X(3872) is a virtual state in the D^D*^ channel, the data on the new enhancement 

Tij- ■ comply with those on the X(3872). In our fits, the mass distribution in the D^D*^ 

o , 

, mode exhibits a peak at 2 3 MeV above the D^D*^ threshold, with a distinctive 



non-Breit-Wigner shape. 



X 

H PACS numbers: 14.40. Gx, 13.25.Gv, 12.39.Mk, 12.39.Pn 



I. INTRODUCTION 

The X(3872) state, discovered by Belle jl| in the B-meson decay, remains the most 
prominent member of the family of "homeless" charmonia, that is those mesons which 
definitely contain a cc pair but do not fit the standard charmonium assignment. The state 
was confirmed then by CDF j^, DO [sl, and BaBar [4]. The charmonium option for the 
X(3872) looks implausible as the state lies too high to be a ID charmonium, and too low to 
be a 2P one j^. This could, in principle, mean that we simply do not understand the spectra 
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of higher charmonia. Indeed, most of the quark model predictions consider charmonia as cc 
states in the quark potential model, with the potential parameters found from the description 
of lower charmonia, with uncertainties coming from proper treatment of relativistic effects. 
Another source of uncertainty is the role of open charm thresholds, the problem which is 
far from being resolved, though the attempts in this direction can be found in the literature 
— see, for example, Refs. 6|, |7|]. In any case, it looks premature to reject the cc assignment 
for the X(3872) on basis of the mass only. However, the further development has revealed 
more surprises. 

The discovery mode of the X(3872) is tt~ J / ip . The observation of the X(3872) in 
the and -n^ -n^ ti^ J / ip {ujJ/ip) modes ^] implies that the X has positive C-parity, and 

the dipion in the t^^-k^ J/ip mode is C-odd, that is it originates from the p. Coexistence of 
the pJ/ip and uJ/ip modes points to a considerable isospin violation. Studies of the dipion 
mass spectrum in X(3872) -k^-k^ J/ip decay establish that only the 1++ or 2"^ quantum 
number assignments are compatible with the data, while all other hypotheses are excluded 
by more than 3(T j^. 

Both 1++ or 2 quantum numbers options for the X(3872) require drastic revisions 
of naive quark potential models, and no alternative explanation of the state in this 
mass region was suggested. On the other hand, it was pointed out in Refs. 
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Hi that 



the DD* system with 1"^+ quantum numbers can be bound by pion exchange, forming a 
mesonic molecule (see also Ref . |l2^ ) ^ . As confirmed by actual calculations [l^ , large isospin 
mixing due to about 8 MeV difference between the D^D*^ and D^D*^ thresholds can be 
generated in the molecular model in quite a natural way. This model was supplied, in 



Ref. 



14| . by quark-exchange kernels responsible for the transitions DD* pj/ip,ujj/ip, 



predicting the uJ/ip decay mode of the X(3872). Note, however, that one-pion-exchange 
as a binding mechanism in the DD* system should be taken with caution, as, in contrast 
to A^A^ case, here the pion can be on-shell, as pointed in ISj, where the ability to provide 
strong enough binding with one-pion exchange was questioned. For the most recent work 



on the implications of the nearby pion threshold see Refs. 



X as quark state we refer to Ref. 



1^ 



16|, 



l3- 



For recent work for the 



and references therein. 



The molecular model has received additional support with the new data on the mass of 



^An obvious shorthand notation is used here and in what follows: DD* = -^{DD* + DD* 



the D meson [19| which yield a very weak binding, 

Mx - M{D^D*^) = -0.6 ± 0.6 MeV. (1) 

I — I 

In the meantime, the Belle Collaboration has reported the first observation [20] of the 
near-threshold enhancement in the D^D^tt^ mode in the decay B KD^D^n^, with the 
branching fraction 

Br{B KD'^D^^) = (1.22 ± O.Sll":^^) • 10^^ (2) 
The peak mass of the enhancement is measured to be 

Mpeak = 3875.2 ± 0.7+?;^ ± 0.8 MeV. (3) 
Obviously it is tempting to relate this new state to the X(3872). However, the average 



value of the X(3872) mass is |2l|] 



Mx = 3871.2 ± 0.5 MeV. (4) 

The central value ([3]) of the D^D^n^ peak mass enhancement is about 4 MeV higher than 
that, which obviously challenges attempts to relate this new state to the X(3872). 

Quite recently, the indication appeared that the Belle result [2^ is likely to be con- 
firmed. Namely, the BaBar Collaboration has reported the preliminary data [22] on the 
B KD^D*^ decay, where the enhancement with the mass of 

M = 3875.6 ± 0.71};^ ± 0.8 MeV, (5) 

was found, in a very good agreement with ([3]). BaBar observes the enhancement in the 
D'^D^TT^ and in the D^D^'j modes, which strongly supports the presence of the D^D*^ 
intermediate state in the decay of the new X. 

If the new BaBar data persist, and the enhancement at 3875 MeV is indeed seen in two 
independent experiments, the possibility should be considered seriously of the presence of 
two charmonium-like states, X(3872) and X(3875), surprisingly close to each other and to 
the threshold. 

However, there exists another, less exotic possibility. Namely, if the X(3872) is indeed 
strongly coupled to the D^D*^ channel, and indeed has l"*"^ quantum numbers, one could 
expect the existence of a near-threshold peak in the D^D*^ mass distribution. In the present 
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paper we perform a phenomenological Flatte-like analysis of the data on the decay B — >■ 
KD^D^Tc^ in the near-threshold region under the assumption of the X —>■ D^D*^ —>■ D^D^-k^ 
decay chain and 1"^"*" quantum numbers for the X. The data on the Ti^Ti^J/ip decay modes 
of the X(3872) are analyzed in the same framework, in order to investigate whether these 
data can accommodate the X(3875) state as a manifestation of the X(3872). 



II. FLATTE PARAMETRIZATION 



In this Section we introduce the Flatte-like parametrization of the near-threshold ob- 
servables. The relevant mass range is between the thresholds for the neutral and charged 
D-mesons. A natural generalization of the standard Platte parametrization for the near- 



threshold resonance 



with 



231] of the D^D*^ scattering amplitude reads 
giKi g2K2 , ,r{E) 



(6) 
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(7) 



and 



5 = M{D+D*-) - M{D°D*°) = 7.6 MeV, 

h = V^inE, = ^-2fiiE, k2 = ^2fi2{E-6), K2 = ^2fi2{6~E). 

Here /ii and fi2 are the reduced masses in the D^D*^ and D^D*^ channels, respectively, and 
the energy E is defined relative to the D^D*^ threshold. In what follows we assume isospin 
conservation for the coupling constants, Qi = 92 = 9- 

The term ir/2 in Eq. (JTj) accounts for non-DD* modes. The X(3872) was observed in 
the n^-K^ J/ip^ 7T^7i^7i^J/ip, and "yJ/ip niodes, with 



Br{X 



H H TT 



Br{X TT+TT- J/^) 

Br{X -iJji)) 



Br[X 



1.0 ±0.4 ±0.3, 



0.14 ±0.05, 



(8) 
(9) 
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reported in Ref. [8j. Thus we assume that r{E) in Eq. ([7j) is saturated by the n~^7i~J/ip and 
71^71^71^ J/ modes and, in accordance with findings of Ref. 8j, the dipion in the tt^tt~ J/ij) 
mode comes from the p whereas the tripion in the n^n^Ti^J/tp mode comes from the u. The 
"jJ/ip channel is neglected due to its small branching fraction ([9]). 

The nominal thresholds for both pJ/tp and ujJ/ip (3872 MeV and 3879 MeV, respectively) 
are close to the mass range under consideration, but both the u meson and, especially, the 
p meson have finite widths, which are large in the scale under consideration. Thus r{E) is 
calculated as 

r{E) = r^+T,-j/^{E) + r^+-^-^oj/^{E), (lo) 
r..,-.„AE)-Ul^^ ^(,„_„^). + r./4. (11) 

with fp and f^j being effective couplings and 



{AP - (m + mj/^y){M^ - (m - mj/^y) 

<iM = V UP (^^^ 

being the cm. dipion/tripion momentum (M = E + M{D^D*^)). 

Now we are in a position to write down the differential rates in the Flatte approximation. 
These are 

dBr{B ^ KD^D*"^) _ 1 gh , . 

dE 2Ti\D{E)\^' ^ ' 

dBr{B ^ Kn+n-J/^) _ 1 E^+^-j/^E) 



dE 2tt \D{E)\ 



and 



dBrjB ^ Kn+Ti-n^J/tP) _ 1 E^+^-^oj/^E) 

dE 2n \D{E)\^ ■ ^ ^ 

We assume the short-ranged dynamics of the weak B ^ K transition to be absorbed into 
the coefficient B. Obviously, the rate (fT^ is defined for > only, while the rates (fTSll and 
f fT6|) are defined both above and below the D^D*^ threshold. 

The formulae fll4l) - fll6l) are valid in the zero-width approximation for the D*-mesons. In 
principle, one could include the finite width of the D*-mesons either analogous to Eqs. fllip 
and ( JT2I) or in a more sophisticated way, as there are interference effects possible in the final 
state, as described in Ref. |2J]. However, the widths of the D* mesons are small. Indeed, 



the total width of the D*^-meson is measured to be 96 ± 22 keV 



2l| . There are no data 
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on the D*" width, but one can estimate the D^vr*^ width of the D*'^ from the data [2l|] on 
charged D*^, which gives r[D*^ D^n^) = 42 keV. The branching fractions of D*^ are 
known 211 ]: 



5r(D*° ^ DV) = (61.9 ± 2.9)%, (17) 

Br{D*^ D^-f) = (38.1 ± 2.9)%, (18) 

so the total D*'^ width can be estimated to be only about 68 keV. The effect of such a small 
width was checked to be negligible in our studies, and we assume the D^D^n^ differential 



rate to be 



dBr(B^KDW)^ ^^^^ 

dE 27i\D{E)\^' ^ ' 

where the branching fraction ( |T7I) is taken into account. 
Analogously we have for the D^D^'-f differential rate 

dE 2n\D{E)\^ ^ ' 

Expressions (1191) and fl20|) neglect final-state interactions; in particular, no DD resonance 
within a few MeV above D^D^ threshold is assumed to exist, and vr-rescattering is neglected. 
The latter is expected to be weak, as a consequence of chiral symmetry [17]. 

III. FLATTE ANALYSIS: PROCEDURE AND RESULTS 

Let us first specify the data used in our analysis. For the tt+tt^J/-?/' mode we use the 



data from the 5- meson decay. These are the ones reported by the Belle [ll] and BaBar [25 ] 
Collaborations. The X(3872) is seen by Belle in the charged 5-meson decay, with 35.7 ±6.8 
signal events, and with the branching fraction [1|] 

Br{B+ K+X)Br{X tt+tt- J/^) = (13.0 ± 2.9 ± 0.7) ■ 10"^ (21) 

n 

The BaBar Collaboration [25'] has observed the X(3872) both in the charged and neutral 
i?-meson decays, with 61.2 ± 15.3 signal events for the charged mode, and only 8.3 ± 4.5 
signal events for the neutral one. The branching fraction for the charged mode was found 
to be 

Br{B- K-X)Br{X vr+vr- J/^) = (10.1 ± 2.5 ± 1.0) ■ 10"^ (22) 
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while the result for the neutral mode is much less certain: a 90% confidence interval was 
estabhshed as 

1.34 ■ 10"^ < Sr(S° ^ K°X)Br{X vt+tt' J/^) < 10.3 ■ 10"^ (23) 

Due to large errors and much smaller number of events, the X(3872) peak in the neutral 

mode looks much less convincing than the peak in the charged mode. 

- □ 
A similar situation takes place for the D^D^tx^ final state. The Belle data [20|] include 

both B+ K+D^D^t:^ and 5° K^D^D^tx^ decays. There are 17.4 ±5.2 signal events in 

the charged mode, with the branching fraction 

Br{B+ K+D'^D\'^) = (1.02 ± 0.31^°;^^) ■ 10-^ (24) 

and 6.5 ± 2.6 signal events in the neutral mode, with 

Br{B^ ^ K^'D^D^^) = (1.66 ± 0.70^°J?) ■ 10^^ (25) 



Data on the B~^ and B^ decays separately are presented in Ref. 26|]. The D^D^ir^ enhance- 
ment appears to be clearly seen in the data on charged B decays while, again, the neutral 
mode displays, within the errors, a much less pronounced peak. 

We conclude therefore that the data on charged and neutral B decays should be analyzed 
separately. The present analysis is performed for the charged mode only. Namely, with the 
Flatte formalism, we attempt to describe simultaneously the TT^n^J/tp mass spectrum from 



the charged mode and the D^D^tt^ spectrum from the B~^ mode, taken from Ref. 261]. 

The branching fractions (l2Tll and (l22l) differ but, within the errors, are consistent with 
each other. In both sets of data, the fitted width of the signal is consistent with the resolu- 
tion, so only the upper limits on the X(3872) width were established: 

rtot(Belle) < 2.3 MeV (26) 

and 

rtot(BaBar) < 4.1 MeV, (27) 

for the Belle and BaBar data, respectively. In view of this discrepancy we prefer to present 
two sets of fits, based on the two aforementioned sets of the Ti^-K^J/ip data. 

The Ti^Ti^ J/ip data are fitted in the interval —20 < < 20 MeV (as before, E is the 
energy relative to the D^D*^ threshold), after subtraction of the full background found in 
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the corresponding analysis. The free parameters of the fit are the short-range factor B and 
the Flatte parameters Ej, g, and fp. The parameter is constrained, in accordance with 
Eq. (IH]) through the condition 

1, (28) 



Rpj/ip 
R 



where 



wj/il) 



20MeV 



20MeV 



dE 



The hmits of integration in Eqs. (129|) and (130|1 are somehow arbitrary but, as most of the 
support of the distributions f|T5|) and f|T6l) comes from within a few MeV around the D^D*^ 
threshold, the uncertainty introduced by the limits of integration is much less than the 
experimental errors in Eq. ([8]). 

The D°Z)°7r° data are fitted in the energy region < E < 20 MeV. Equation (|T9ll 
describes the production of the D^D^tx^ mode via the X-resonance, while the DD* pairs are 
known to be copiously produced in the B ^ K decay in a non-resonant way. Besides, the 
D^D^-K^ final state could come from non-D°Z)*° modes like, for example, B K*D^D^. 
Therefore, we are to make assumptions on the background. 

The background in Refs. j2o| and 26 1 is mostly combinatorial, and this part, given explic- 
itly in the publications, was subtracted prior to the analysis. For the rest of the background 
it is not possible to separate the contributions of the D^D*^ and the D^D^tt^ due to a limited 

n 

phase space [20|. So we work under two extreme assumptions for the background. In Case 
A we consider the D^D^n^ background as unrelated to the D^D*^ channel, while in Case B 
we assume that all the D^D^n^ events come from the D^D*^ mode. The background was 
evaluated by fitting the Belle data off-peak (25 < E < 50 MeV). In Case A the background 
function is assumed to be proportional to the three-body D^D^tt^ phase space R3 oc E]~,j-,^, 
where -Eddtt = E + m£)*o — m^io — m^o. Then the total B KD^D^n^ differential rate is 

dBr^iB ^ KD^D^Ti^) B gh 2 

with ca as fitting constant. In Case B the background function is proportional to the 
two-body D^D*^ phase space R2 oc ki (see the definition below Eq. ([7])). Then the signal- 
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TABLE I: The set of the Flatte parameters for the best fits to the Behe data Ref. II] and 20l |. 



Fit 


9 


fp 




Ef, MeV 







^Belle 


0.3 


0.0070 


0.036 


-11.0 


11.0 






0.3 


0.0086 


0.046 


-10.9 


8.9 


-144° 



background interference is to be taken into account: 



(32) 



with the relative phase (p and cb being fitting constants. 

The differential rates are translated into number-of-events distributions as follows. There 
are about 36 signal events in the Belle data, which corresponds to the branching fraction of 
about 1.3 ■ 10~^ (see Eq. (I2T!) ). Thus the number-of-events per 5 MeV distribution for the 
Ti^-K^ J/ip mode is given by 

N'^i*(E) = 5 IMeV] (^^) "^'''^ K.^-' J hi') , (33) 

For the BaBar data, with 61 events and the branching fraction of about 1.02 ■ 10^^ (see 
Eq. ([22])), we have 

Nl'^:(E) = 5 IMeV] {j^) iBr(B ^ K.-.- J 1^) 

As to the D^D^n^ mode, the Belle Collaboration states to have 17.4 signal events in the 
oha.,ed „„de H. wh.cK con.e.po.d. to the b..a.ch,.g f.act.ou » of about 1.02 • 10-. 

The number-of-events distributions per 4.25 MeV for the D^D^tt^ mode is calculated as 



- (E) = 4,25|MeVl ^^^^ j ^— (35) 

The best fit to the TT^TT^J/tp data alone requires a vanishing value of the DD* coupling 
constant, g = 0, so that such solution cannot accommodate the D^D^tx^ enhancement as a 
related phenomenon. To describe both -n^-n^ J and D^D*^ modes we are to compromise 
on the TT^TX~ J/ip line-shape. 

It appears that a decent combined fit can be achieved only if the -k^-k^ J/ip distribution 
is peaked exactly at the D^D*^ threshold, with the peak width (defined as the width at the 
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TABLE II: The set of the Flatte parameters for the best fits to the BaBar data of Ref. 25| and 



the Behe data of Ref. 



2C| 



Fit 


9 






Ef, MeV 


10^ 




^BaBar 


0.3 


0.0042 


0.021 


-8.8 


11.4 




-^BaBar 


0.3 


0.0056 


0.027 


-8.8 


8.9 


-153° 



peak half-height) close to the upper limits given by Eq. (|26|) or fl27|) . The values of the 
coupling g were found to be of the order of magnitude or larger than 0.3. Finally, the fits 
exhibit the scaling behaviour: they remain stable under the transformation 

g^Xg, Ef^XEf, f,^\f„ ^ \f^, B ^ \B, (36) 

with tiny variations of the phase in the Case B. 

In Tables [H [III we present the sets of the best fitting parameters — for both Case A 
and Case B and for g = 0.3 — for the Belle (Table [T]) and BaBar (Table [TTl) data on the 
71^71^ J/ip mode and for the Belle data for the D^D^tt^ mode. To assess the quality of the 
fits we calculate the n^n^J/ip distributions integrated over the 5 MeV bins, as in Refs. [1| 
and |25], and the D^D^tt^ distributions integrated over the 4.25 MeV bins, as in Refs. 20 1 
and [26]. The results are shown at Fig. [1] together with the experimental data. 

The above-mentioned scaling behaviour does not allow one to perform a proper fit with 
the estimate of uncertainties in the parameters found. Indeed, the parameters of the best 
fits found for the values of coupling constant g larger than 0.3 differ only by a few % from 
the ones given by the scaling transformation (J36l) . and the corresponding distributions are 
very similar to those given at Fig. [H 

As seen from the figures, acceptable fits require the D^D*^ differential rate to be peaked 
at around 2 -h 3 MeV above the D^D*^ threshold. The scattering length in the D^D*^ 
channel which follows from the expression ([6]) of the D^D*^ scattering amplitude, is given 
by the expression 



V2iI^ + 2Ef/g + zr{0)/g ^^^^ 



iVW + 2Ef/g)^ + r{0)yg^' 
and is calculated to be 



(-3.98 - iO.46) fm. Case Aeeiie 

(38) 

(—3.95 — ^0.55) fm. Case Beeiie, 
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FIG. 1: Upper plots: Our fits to the differential rates for tlie tt+tt J/tp channel measured by 
Belle [ll and BaBar 25] using prescription A and B (see Eqs. ([3T]) and ([32]) 1. Lower plots: Cor- 



2g]. The 



responding fits for the differential rates in the D^D^ir^ channel measured by Belle 
distributions integrated over the bins are shown in each panel as filled dots, experimental data as 
filled squares with error bars. 



(39) 



and 

(—3.10 — z0.16) fm, Case AsaBar 
(-3.10 - iO.22) fm, Case BBaBar- 
The real part of the scattering length for all the fits appears to be large and negative, and the 
imaginary part is much smaller. This, together with the beautiful cusp in the -k^-k^ J /ip mass 
distribution, signals the presence of a virtual state in the D^D*^ channel. The cusp scenario 
for the -n^Ti^J/ij) excitation curve in the X(3872) mass range was advocated in Ref. ^!\ - 



The X 
model 



3872) as a virtual DD* state was found in the coupled-channel microscopic quark 
1 



A large scattering length explains naturally the scaling behaviour of the Flatte param- 
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eters. Such kind of scaling was described in Ref. [28| in the context of hght scalar mesons 
properties: the scaling behaviour occurs if the scattering length approximation is operative. 
In the case of X the situation is more complicated, as there are two near-threshold channels, 
neutral and charged. Nevertheless, if it is possible to neglect the energy E in the expression 
for the Flatte denominator D{E) then, as seen from the expression ([6]), scaling for the 
scattering amplitude indeed takes place. If the factor B obeys the scaling transfor- 
mation, the differential rates f|T^ - f|T6|) also exhibit the scaling behaviour. Note that, if the 
energy dependence of the charged D^D*^ and non-DD* channel contributions is neglected 
as well, this corresponds to the scattering length approximation, and neglect of the effective 
radius term. 



IV. DISCUSSION 

Our analysis shows that the large branching fraction ([2]) implies the X to be a virtual 
£)0£)*o ^loXe, and not a bound state. We illustrate this point by calculating the rates (fT^ 
and (fTSl) for the set of the Flatte parameters (fit C) 

g = 0.3, Ef = -25.9 MeV, fp = 0.007, = 0.036, B = 1.32 ■ 10"^ (40) 

The values of the coupling constants coincide with those of the fit Aeeiie, while the parameter 
Ef is chosen to yield the real part of the scattering length to be equal in magnitude to the 
one evaluated for the given fit ^Beiie, but positive: d = (+3.98 — ■jO.46) fm. The parameter 
B for this set yields the same value of the total branching fraction for the -k^-k^ J /ip mode as 
the fit A-QeWe- The t^^-k^ J /iIj and D^D^tt^ rates are shown in Fig. [2l together with the rates 
obtained for the case ^Beiie (without background). The new curve (dashed line in Fig. [2]) 
displays a very narrow peak in the ti^-k^ J /ip distribution, corresponding to the D^D*^ bound 
state, with binding energy of about 1 MeV (there is no corresponding peak in the D^D^tt'^ 
distribution as the finite width of the D*^ is not taken into account in our analysis). Note 
that the h^tt' J /if) rates (Fig. [2]) are normalized to give the branching ratio 1.3 ■ 10~^, which 
requires the coefficient B to be much larger for the virtual state than for the bound state. 
As a result, the D^D*^ rate is much smaller for the bound state, as seen from Fig. O 
Obviously, the difference between the bound-state and virtual-state cases for the ratio 
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FIG. 2: The differential rates for the ir'^Tr J/t/j (first plot) and D^D*^ (second plot) for the fits 
^Beiie (solid curves) and C (dashed curves). 

is driven by the strength of the bound-state peak, as discussed in Ref. [3], where the 
scattering length approximation was used to describe the X(3872). Following Ref. 29||, let 
us write down the scattering length in the D^D*^ channel as 

1 



- ■• (42) 

Then, in the scattering length approximation, the 7T~^7r~J/ip differential rate is proportional 
to the factor 



Ire + {ki+limf 



(7r 



E > 0, 



E < 0, 



(43) 



while the D^D*^ rate is proportional to 



l're + ih+l^n.r ^^^^ 

The line-shape for the D^D*^ channel does not depend on the sign of ''jre- The same is 
true for the -k^ti^J/iI) line-shape above the D^D*^ threshold while, below the threshold, 
the line-shapes differ drastically: in the bound-state case there is a narrow peak below 
threshold, and in the virtual-state case a threshold cusp appears. 
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For > and the expression (H3i) becomes a (5-function (see Ref. |29[|): 

-7.e5(i? + 7'e/(2/il)). (45) 

Then the total rate does not depend on 7jm, if it is small enough. This simply means that, 
for 7jm = 0, we have a real bound state, which is not coupled to inelastic channels. In 
contrast to the bound-state case, for the virtual state, the rate (H3|) tends to zero with 
7im 0, while the D^D*^ rate does not vanish in such a limit. So it is possible, adjusting 
7im, to obtain large values of the ratio (HB . 

Exactly the same situation is encountered in our fit: we need g > 0.3 for the fit to 
be reasonable and, in this scaling regime, as soon as we have a positive real part of the 
scattering length, the ratio (jHj) becomes small while, with a negative real part, we get 
a solution compatible with the data. The large branching fraction ([2]) was identified in 
Ref. |30|] as a disaster for the molecular model of the X(3872). Indeed, the bound-state 
molecule decay into D°Z)°7r° is driven by the process D*^ D^tt^ which gives the width 
of order 2T(D*^ D'^n'^) (up to the interference effects calculated in Ref. |2j] which, for 
the bound-state case, cannot be neglected anymore and should be taken into account). The 
main decay mode of the X is n^n^J/ip because the phase space available is large. This is 
confirmed by model calculations of Ref. [mI yielding 

Br(X ^ D^D^TT^) 

7^ ~ 0.08, 46 

5r(X ^ vr+vr- J/^) ' ^ ^ 

in a strong contradiction with data. 

The estimate P6|) describes the decay of an isolated bound state. However, the suppres- 
sion is more moderate as, in B-decay, the continuum contribution is also to be considered. 
The bound-state contribution would be zero in the zero-width approximation for D*^, while 
the D^D*^ continuum contribution remains finite if the D*^ width is neglected. However, if 
the X is a bound state, the continuum contribution is not large (see Fig. [2]), 

Br(X D^D^Tt^) 

Such a small rate would remain unnoticed against the background. So, in practice, the 
bound-state X(3872) would reveal itself only as a narrow peak below threshold, with a very 
small rate (see Eq. (l46l) ). In contrast to this we get for the virtual state 

Br{X D^D^n^) 

Br{X ^ 71+71- J/^) ^ ^ 
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In our analysis, the X appears to be a virtual state in the D^D*^ channel. This does not 
contradict the assumption Qi = g2 = 9 employed in the analysis. The latter means that the 
underlying strong interaction conserves isospin, and all the isospin violation comes from the 
mass difference between charged and neutral DD* thresholds. No charged partners of the 
X are observed, so it is reasonable to assume that the strong attractive interaction takes 
place in the isosinglet DD* channel. 

We do not specify the nature of this attractive force. It is known that in the one-pion- 
exchange model for the X, the force is attractive in the isosinglet channel, and is repulsive 
in the isotriplet one. However, as was already mentioned, the doubts were cast in 15|] on the 
role of one-pion-exchange in the DD* binding, and it was advocated there that the X may 
fit the 2^ Pi charmonium assignment if the coupling to DD* channel is taken into account. 
In such a scenario the strong binding force obviously takes place in the isosinglet channel. 

We note, however, that, with the Flatte parameters found, one can make a definite 
statement: whatever the nature of the X(3872) is, the admixture of a compact cc state in 
its wavefunction is small. Both large scattering length and the scaling behaviour of the DD* 
amplitude are consequences of the large value of the coupling constant of the state to the 



DD* channel. As shown in Ref. 



3l|, this points to a large DD* component and a dynamical 



origin of the X. Although formulated for quasi-bound states in Ref. |31|] the argument can 
also be generalized to virtual states. To clarify the connection between effective coupling 
and the nature of the state observe that the two-point function g{s) for the resonance can 
be written as 

where M is the physical mass of the resonance and = S(s)— ReS(M^) is the self-energy 
responsible for the dressing through the mesonic channels. In the near-threshold region the 
momenta involved are much smaller than the inverse of the range of forces. As a result one 
may neglect the s-dependence of the real part of S and replace its imaginary part by the 
leading terms 

where the sum is over near-threshold channels, and the contributions of distant thresholds 
are absorbed into the renormalised mass M. Nonrelativistic reduction of Eq. ( l50l) immedi- 
ately yields the Flatte formula ([6]). Thus the Flatte parameter Ef acquires clear physical 
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meaning: the quantity M{D^D*'^) + Ef is the physical mass of the resonance, renormahsed 
by the couphng to the decay channels. 

Now, if the couplings Qi are small, the distribution for the resonance takes a standard 
Breit-Wigner form, and the scattering length is small. Correspondingly, the state is mostly 
cc, with a small admixture of the DD* component. If the couplings are large, the terms 
proportional to Qiki control the denominator in Eq. (l50|l . the Breit-Wigner shape is severely 
distorted, the scattering length approximation is operative, and the mesonic component 
dominates the near-threshold wavefunction. 

Formulated differently: if the couplings are large, the properties of the resonance are 
given mainly by the continuum contribution — which is equivalent to saying it is mostly 
of molecular (dynamical) nature. It should be stressed that this kind of reasoning can only 
be used, if the resonance mass is very close to a threshold, for then the contribution of 
the continuum state is dominated by the unitarity cut piece which is unique and model 
independent. This argument is put into more quantitative terms in Ref. [3l|]- It is also 
important to note that our analysis does not allow for any conclusion on the mechanism that 
leads to the molecular structure. On the level of the phenomenological parametrisations used 
here a molecule formation due to t-channel exchanges and due to short-ranged s-channel 
forces {cc-DD* mixing) would necessarily lead to the same properties of the state, once the 
parameters are adjusted to the data. 



V. SUMMARY 

n 

In this paper we present a Flatte analysis of the Belle data [20] on the near-threshold 
enhancement in the D^D^tc^ mode. We constrain the Flatte parametrization with the data 
on the X(3872) seen in the n^Tr^J/tp and tt+tt^tt^J/?/' modes. With such constraints the 
new state can be understood as a manifestation of the well-established X(3872) resonance. 

We showed that the structure at 3875 MeV can only be related to the X(3872), if we 
assume the X to be of a dynamical origin, however, not as a bound state but as a virtual 
state. The situation is then similar to that of nucleon-nucleon scattering in the spin- 
singlet channel near threshold: in contrast to the spin-triplet channel, where there exists 
the deuteron as a bound state, the huge scattering length in the spin-singlet channel — 
about 20 fm — comes from a near-threshold virtual state. The attractive interaction is just 
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not strong enough to form a bound state in this channel as well. 

The line-shape in the D^D*^ mode appears to differ substantially from the one extracted 
previously from the Belle data directly. It peaks much closer to the D^D*^ threshold, though 
the overall description of the data looks quite reasonable within the experimental errors. 

It is the Tx^-K^J/ip line-shape which, in our solutions, differ drastically from the one 
described by a simple Breit-Wigner form. We found a threshold cusp, with a width close to 
the limits imposed by the data analysis. While the data currently available allows for such 
a line-shape, a considerable improvement in the experimental resolution could confirm or 
rule out this possibility. In the meantime, we urge to perform an analysis of the data on the 
D^D^Tt^ final state with Flatte formulae given in Eqs. (fT^ - (fT6i) . 

Equally important is the Flatte analysis of the D^D^'-y data j^]: if the structure in the 
D^D^TT^ is indeed due to D^D*^ and is indeed related to the X(3872) as a virtual state, one 
should observe an enhancement in £)°Z)°7 similar to the one seen in the D^D^tt^. The phase 
space available in this final state is larger than that in D^D^n^, so it is easier to separate 
the contributions of D^D*^ and D^D^'-f to the peak. The D^D^j enhancement would be 
described with the Flatte formula ( l20l) and, up to background and possible FSI effects, the 
ratio of branching fractions would be 

Br{X -> D^D^n^) 

Br{X D^D^^) ^-^^ ^^^> 

The most interesting situation would happen if, due to an improved resolution in the 
Ti^-K^ J/ip mode, the combined Flatte analysis of the -k^-k^ J/ip, D^D^tt^, and D^D^'-y data 
fails to deliver a self-consistent result. Such a situation would point to the new X(3875) 
state being completely unrelated to the X(3872). 
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